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^D Abstract 

(N 

This paper is concerned with the study of a model case of first order Hamilton-Jacobi 
^ equations posed on a "junction", that is to say the union of a finite number of half-lines with 

. ^ a unique common point. The main result is a comparison principle. We also prove existence 

1^1 and stability of solutions. The two challenging difficulties are the singular geometry of the 

f-y-j domain and the discontinuity of the Hamiltonian. As far as discontinuous Hamiltonians are 

^vj concerned, these results seem to be new. They are applied to the study of some models arising 

in traffic flows. The techniques developed in the present article provide new powerful tools 
for the analysis of such problems. 
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^ 1 Introduction 

> 



In this paper we are interested in Hamilton-Jacobi (HJ) equations posed on a one dimensional 
^o domain containing one single singularity. This is a special case of a more general setting where HJ 

CO equations are posed in domains that are unions of submanifolds whose dimensions are different 

p I [5]. An intermediate setting is the study of HJ equations on networks, see in particular [T|. We 

f^ will restrict ourselves to a very simple network: the union of a finite numbers of half-lines of the 

^-H plane with a single common point. Such a domain is referred to as a junction and the common 

'"^ point is called the junction point. We point out that getting a comparison principle is the most 

^ difficult part in such a study; it is obtained in [1 for similar special networks (bounded star-shaped 

ones). Our motivation comes from traffic flows. For this reason, it is natural to impose different 

^ . dynamics on each branch of the junction. Consequently, the resulting Hamiltonian is by nature 

C^ discontinuous at the junction point. Together with the singularity of the domain, this is the major 

technical difficulty to overcome. The analysis relies on the complete study of some minimal action 

(or metric) related to the optimal control interpretation of the equation [33|, I16j . We prove in 

particular that this minimal action is semi-concave by computing it. 

We first present the problem and the main results in details. Then we recall existing results 
and compare them with ours. 
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1.1 Setting of the problem 

In this subsection, the analytical problem is introduced in details. Wc first define the junction, 
then the space of functions on the junction and finally the Hamilton-Jacobi equation. 

The junction. Let us consider A^ > 1 different unit vectors Ci G R'^ ior i — 1, ...,N. We define 
the branches 

J, = [0, +c») • e„ j; = J,;\{0}, z = l,...,iV 

and the junction (see Figure nl) 

j= u A- 

i=l,...,N 

The origin a:: = is called the junction point. It is also useful to write J* = J \ {0}. For a time 




Figure 1: A junction 

T > 0, we also define 

Jj, = (0, T) X J. 

The reader can remark that we chose to embed the junction in a two-dimensional Euclidean space. 
But we could also have considered an abstract junction, or we could have embedded it for instance 
in a higher dimensional Euclidean space. We made such a choice for the sake of clarity. 

Space of functions. For a function u: Jj- — >■ M, we denote by u' the restriction of u to (0, T)xJi. 
Then we define the natural space of functions on the junction 

Cl{JT)^{ueC{JT), u' e C\{0,T) X Ji) for i = l,...,N}. 

In particular for u e C'KJt) and x ~ XiCi with Xi > 0, we define 

ut{t,x) ^ —{t,Xiei) and ul{t,x) ^ —{t^XiCi). 

Then we set 

Ux{t,x) — ul.{t,x) if X ^ 0, 

u,(t,0) = K(i,0))j=i,...^w if x = 0. 
HJ equation on the junction. We are interested in continuous functions u: [0,r) x J ^- M 



which are viscosity solutions (see Definition 3.2) on Jy of 



Ut + H{x,u,)^0 (1.1) 

submitted to an initial condition 

u{0,x) ^uq{x), X e J. (1.2) 



Because of the optimal control problem we have in mind (see Subsection |1.1| below) , we restrict 
ourselves to the simplest case of discontinuous Hamiltonians; precisely, we consider 



H,{p) for peR a xeJ* 

niaxj=i,...^jv^r(Pi) for P ^ {Pi, ■■■^Pn) '^R^ if x^O 



where Hi are convex functions whose Legendre-Fenchel transform is denoted Li. We recall that 

Hiip) = L*{p) = sup {pq - L,{q)) 

and 

H.=L*. (1.3) 

We also consider 

H-ip,)^snpip,q-LM) (1-4) 



Therefore equation (1.1) can be rewritten as follows 



ul + Hiiul)^0 on (0,r)xJ* for i^l,...,N, 

Uf +max,;=i,...,7vi?~«) = on (0,r) x {0}. 



(1.5) 



The optimal control framework. In this paragraph, we give an optimal control interpretation 
PSI [71 2] of the Hamilton-Jacobi equation. We define the set of admissible controls at a point 
X £ J hy 



^^""^"1 U=i,...,ArK+e, if x = 

For {s,y), {t,x) G [0,r] x J with s < t (the case s — t being trivial and forcing y = x), we define 
the set of admissible trajectories from (s, y) to {t, x) by 



Ais,y;t,x) = {X e W'''{[s,t];R''): 



X{t)£J for all t £ {s,t), 

X{t)£A{X{t)) fora.e. T£{s,t) ). (1.6) 

xls) = y and X{t) = x 



For P = pci G A{x) with p G M, we define the Lagrangian on the junction 



^(-•^)HtS ;[ ^=o' '-"■> 



with 



Lo[p) = . mwL Lj{p). 

j = l.,...,N 

The reader can be surprised by the definition of L{x, P) for x = 0. In fact, if one considers only 
trajectories that do not stay at the junction point, then the ones staying at are approximated 
by those staying very close to on a branch i E Iq and moving "slowly" {X ~ 0). 

1.2 Main results 

We make the following assumptions: 

(AO) The initial data uq is Lipschitz continuous. 



(Al) There exists a constant 7 > 0, and for all i = 1,...,N, there exists C^(IR) functions Li 
satisfying L'- > 7 > 0, such that (|1.3|) and (|1.4|) hold. 



Theorem 1.1 (Existence and uniqueness). Assume (AO)-(Al) and let T > 0. Then there exists a 
unique viscosity solution u of \1.1^ - \1^ on Jx in the sense of Definition 3.2, satisfying for some 
constant Ct > 

\u{t,x) — Uf){x)\ < Ct for all {t,x) £ Jt- 
Moreover the function u is Lipschitz continuous with respect to (t, x) on Jt- 

On one hand, we will see below that the existence of a solution can be obtained with Perron's 



method under weaker assumptions than (Al) (see Theorem A.2|. On the other hand, we are able 



to get uniqueness of the solution only under assumption (Al) and this is a consequence of the 
following result: 

Theorem 1.2 (Comparison principle). Assume (AO)-(Al). Let T > and let u (resp. v) he a 
sub-solution (resp. a super-solution) of [Lw-(L2) on Jt in the sense of Definition 3.2 We also 
assume that there exists a constant Ct > such that for all {t, x) £ Jt 

u{t,x)<CT(l + \x\) (resp. v{t,x) > -Ct{1 + \x\)) . 

Then we have u < v on Jt. 

In order to prove this strong uniqueness result, we will use in an essential way the value function 
associated to the optimal control problem described in Subsection |1.1[ for t > 0, 

nt 



Uoc{t,x) 



inf 

yeJ, XeA(0,y;t,x) 



uoiy) 



LiX{T),X{T))dT 



(1.8) 



where L is defined in (1.7 1 and A{0,y;t,x) is defined in (1.6 1 



Theorem 1.3 (Optimal control representation of the solution). Assume (AO)-(Al) and let T > 0. 
The unique solution given by Theorem \l.l\ is u = Uqc with Uoc given in (uTq). Moreover, we have 
the following Hopf-Lax representation formula 



,{t,x) = inf 1^0(2;) +'L>{0,y;t,x)} 



(1.9) 



with 



2?(0, y; t, x) = inf 

XeA{0,v:t,x) 



L{X{T),X{T))dl 



The comparison principle is obtained by combining 

• a super-optimality principle for counterrevolutions v^ which implies v > Uqc; 

• a direct comparison result with sub solutions w, which gives Uqc > u. 

We finally have the following result which shed light on the role of the junction condition (see the 
second line of (1.5)). 



Theorem 1.4 (Comparison with continuous solutions outside the jmiction point). Assume (AO)- 
(Al) and let T > 0. Let u G C([0,T) x J) be such that u(0, •) = uo and for each i e {l,...,Af}, 
the restriction u' of u to (0,T) x Ji is a classical viscosity solution of 



ul + H,{ui)^0 on (0,r)xJ*. 



Then u is a sub-solution of ['1.1)-['1.2) on Jt in the sense of Definition 3.2 



and u < Ur, 



An immediate consequence of Theorem |1.4| is the fact that the optimal control solution Uqc is 
the maximal continuous function which is a viscosity solution on each open branch. 

We apply in Section [2] our HJ approach to describe traffic flows on a junction. In particular, 
we recover the well-known junction conditions of Lebacque (see |22| ) or equivalently those for the 
Riemann solver at the junction as in the book of Garavello and Piccoli [18]; see also |19) . 



1.3 Comments 

We already mentioned that the main difficulties we have to overcome in order to get our main re- 
sults are on one hand the singular geometry of the domain and on the other hand the discontinuity 
of the Hamiltonian. 



Discontinuity. Several papers in the literature deal with HJB equations with discontinuous 
coefficients; see for instance [3[31llHll3ZllMl[nil3Hl[I3[Il]- Note that in these works the optimal 
trajectories do cross the regions of discontinuities (i.e. the junction in the present paper) only on 
a set of time of measure zero. In the present paper, the optimal trajectories can remain on the 
junction during some time intervals, and the results cited above do not apply then to the problem 
studied here. 

On the other hand, the analysis of scalar conservation laws with discontinuous flux functions 
has been extensively studied, we refer to [35J[21[T^ and references therein. We also point out that 
a uniqueness result is proved in M in the framework corresponding a junction with two branches 
|17j . To the best of our knowledge, in the case of junctions with more than two branches, there 
are no uniqueness result. Moreover, the link between HJB equations and conservation laws with 
discontinuous has been seldom investigated |29j . 

The main differences between the study in [T] and the one carried out in the present paper lie in 
the fact that in |lj the Lagrangian can depend on x and is continuous with respect to this variable, 
while we consider a Lagrangian which is constant in x on each branch but can be discontinuous 
(with respect to x) at the junction. Wc point out that we cannot extend directly our approach to 
Lagrangians depending on x since we use extensively the representation formula "a la Hopf-Lax" . 
In order to generalise results in this direction, the semi-concavity of the "fundamental solution" 
V should be proved without relying on explicit computations. This question is very interesting 
but is out of the scope of the present paper. 

Networks. It is by now well known that the study of traffic flows on networks is an important 
motivation that give rise to several difficulties related to scalar conservation laws with discontinu- 
ous coefficients. This topic has been widely studied by many authors, see for instance [lOl UHl [15] 
and the references therein. 

However, the study of HJB equations on networks has been considered very recently; the reader 
is referred to |31llll| where Eikonal equations are considered. A more general framework was also 



studied in [Tl[2] where a definition of viscosity solutions on networks, similar to Definition 3.2 has 
been introduced. 

It would be interesting to extend the results of the present paper to more general networks 
but the obstacle is the same than the one to be overcome if one wants to deal with Lagrangians 
depending on x: for a general network, the complete study of the fundamental solution is probably 
out of reach. This is the reason why we only consider the very specific case of a junction in order 
to be able to overcome the difficulty of the discontinuity of the Lagrangian. 

The optimal control interpretation. As explained above, the comparison principle is 

proved by using in an essential way the optimal control interpretation of the Hamilton- Jacobi 
equation. The use of representation formulas and/or optimality principles is classical in the study 
of Hamilton- Jacobi equations |57l [3S1 |3S1 [201 HI] • More specifically, it is also known that a "metric" 
interpretation of the Hamilton- Jacobi equation is fruitful [33 . Such an interpretation plays a 
central role in the weak KAM theory [1^ . 

As far as our problem is concerned, we are not able to prove uniqueness of viscosity solutions by 
using the classical techniques of doubling variable. The idea used here is based on the equivalence 
between the viscosity super-solution and the super-optimality principle (also known as weak- 
invariance principle), and by using representation formulas for the viscosity sub solutions. This 



representation seems to be new for HJB equations with discontinuous coefficients, see for instance 

We wouid iilce next to be a bit more precise. Tlie teclmicai core of tire paper lies in Theorem |4.2[ 
This result implies that the function 

'D{s,y;t,x) = {t - s)!)!, 



t-s' t- 

is semi-concave with respect to {t,x) and (s, y) and, if there are at least two branches {N > 2), 

that T> satisfies 

Vt + H{x,V,) =0, 

-V, + H{y,-Vy) =0 



(in a weak sense made precise in the statement of Theorem 4.2 1. In the case where the Lagrangians 
coincide at the junction point {Li{0) = ... = Ljv(O)), it turns out that the restriction T>Q^{y,x) of 
Vq to Jj X Ji belong to C^{Jj x Ji) and is convex. A more general case is considered in this paper: 
Lagrangians can differ at the junction point and in this case, the functions 2?q* are not convex nor 
C^ anymore for some («,j). Let us point out here that the assumptions on the Hamiltonian Hi, 
and in particular the fact that it does not depend on the space variable x, plays a crucial role to 
establish the properties satisfied by the minimal action function V. 

Generalization and open problems. Eventually, we briefly mention natural generalizations 
of our results and some important open problems. First of all, it would be natural to extend 
the results of this paper to more general setting where the Hamiltonians Hi depend on the space 
variable x. Moreover, it would be interesting to consider general networks with several junction 
points. Dealing with non-convex and non-coercive Hamiltonians is quite challenging and would 
require first to have a direct proof of the comparison principle which does not need to go through 
the interpretation of the viscosity solution as the value function of an optimal control problem. 

Organization of the article 

Section [2] is devoted to the application of our results to some traffic flow problems. In particular, 
the HJ equation is derived and the junction condition is interpreted. In Section [3| the definition 
of (viscosity) solutions is made precise. In Section Hi the first important properties of optimal 
trajectories are given. Section [5] is devoted to the proof of the main results of the paper. In 
particular, the comparison principle is proved by proving a super-optimality principle and by 
comparing sub solutions with the solution given by the optimal control interpretation of the 
equation. Section |6] is devoted to the proof of the technical core of the paper, namely the existence 
of test functions for the minimal action associated with the optimal control interpretation. 

Notation 

Distance and coordinates in the junction. We denote by d the geodesic distance defined 
on J by 

, , f |x — j/| if x,y belong to the same branch Ji for some i, 

|_ |x| -f- |y| if x,y do not belong to the same branch. 

For X E J, B{x, r) denotes the (open) ball centered at x of radius r. We also consider balls 
B((t, x), r) centered at (t, a;) G (0, -f cx)) x J of radius r > 0. For x e J, let us define the index i{x) 
of the branch where x lies. Precisely we set: 

Up to reordering the indices, we assume that there exists an index fco e {1, ..., N} such that 

Lo{0) = ii(0) = • • • = ife,(0) < Lfc„+i(0) < • • • < Ln{0). (1.10) 



We also set 

/o-{l,..,fco} and In=^{1,...,N}. 

Functions defined in J'^. For a function ip defined on J x J, we call ip''^ its restriction to Ji x Jj. 
Then we define the space 

Cl{J^) = {ipeC{J^), ^'^eC\j,xJj) for all z,je/^}. 

We also call for x — XiCi with Xi > and y — yjEj with yj > 

d d 

dl^{x,y) = T^v''{x^e,,y) and d^ipix^y) = —f'^ {x,yjej) 



and 



and similarly 



diip{x,y) if xeJ*, 



dMx.y)-[ (ai^(;,y))^^^^ 



N 



if x = 



We also set 



« ^ ^ / 9lipix,y) if 2/e J*, 

5.^(-,y) = | (a^^(x,,))^.^^^^^^ if , = 0. 

xdj;ip{x,y) = Xidlip{x,y) if x e J^, 
ydy(p{x, y) = yjO^ifix, y) if y G J,. 

2 Application to the modeling of traffic flows 

In this section we present the derivation of the Hamilton- Jacobi formulation of traffic on a junction. 
We also discuss the meaning of our junction condition in this framework and relate it to known 
results. 

2.1 Primitive of the densities of cars 

We consider a junction (represented on Figure ^ with m > 1 incoming roads (labeled by the 




Figure 2: A traffic junction 

index i = 1, ..., m) and n > 1 outgoing roads (labeled by j = m + 1, ...,m + n = N). This means 
that the cars move on the incoming roads in the direction of the junction and then have to choose 
to go on one of the n outgoing roads. We assume that the proportion of cars coming from the 
branch i — 1, ...,m is a fixed number 7^ > (which may be not realistic for to > 2), and that the 
proportion of cars going on each branch j e {m + 1, ..., to + n} is also a fixed number 7^ > 0. We 
also assume the obvious relations (for conservation of cars) 

Yl '^' = 1 ^'^^ Yl ^i = 1 

2— l,...,?n j—m-{-l....,m+n 

We denote by p^{t,X) > the car density at time t and at the position X on the branch k. In 
particular, we assume that the traffic is described on each branch fc by a fiux function /'^ : M — >■ M. 
We assume 

each function /'^ is concave and has a unique maximum value a,t p — p'^ > (2.1) 



The typical example of such flux function is given by the LWR model (Lighthill, Whitham [25] 
and Richards [SIT), with 

f{p) = pv{p) with the velocity v{p) = Kiax(l ~ p/Pmax) (2.2) 

where Vmax and /Omax are respectively the maximal velocity and the maximal car density in the 
model). In this model the critical car density pc where / is maximal, is equal to ^Pmax. 
We assume that the car densities are solution of non linear transport equations: 

Pl + irip'))x = 0, X<0, for t^l,...,m 

pi + {P{p^))x=0, X>0, for j^m + l,...,m + n ^'' 

where we assume that the junction point is located at the origin X — 0. 

We do not precise yet the junction condition at X = 0, and we now proceed formally to deduce 
the Hamilton-Jacobi model of such a junction. For a function g to determine, let us consider the 
functions 



W{t,X)=g{t) + ^ J^^ p'{t,Y) dY, X<0, for i = l,...,m. 



7* Jo 



X 



[ W{t,X)=g{t) + ^J^ p^{t,Y)dY, X>0, for j=m + l,...,m + n. 
Then we can compute formally for j = m + 1, ...,m + n 
m ^g'{t) + ^!^p^,{t,Y)dY 

^9'{t)-^J,^ifHp'it,Y)))xdY 



(2.4) 



= -^pip^it,X))+g'it) + ^p(p,(t,0+)) 



This shows that for j = m + 1, ...,m + n 

ui + ^ri^,uj,) = h^t) (2.5) 

where 

Remark that we can show similarly that ( |2.5| is still true for the index j replaced by i = 1, ...,to 
with 

K\t)=g'{t) + -J\p\t,Q-)). 



In particular, this shows (at least when the quantities in (2.5) are well defined) that we can choose 
g such that the total flux —g'{t) is given by 

-g'it)^ J2 fip'it^o-)) (2.6) 

2— l,...,m 

and then we have 

hf'(t) = Q for i = l,...,m + n 

if and only if 

f(p*(t,0-))=y(-g'(t)) for z-l,...,m 

P(pJ(i,0+)) =7^(-g'(i)) for j = m + 1, ...,m + n ^ ' 

which is exactly the expected condition which says that the proportion of incoming cars going in 
the junction from the branch i is 7* and the proportion of cars getting out of the junction which 
choose to go on the branch j is 7-' . 



Let us notice that if we choose the initial condition 17(0) = 0, then we deduce from (2.6) that 
we have with I — 1, ...,m + n 



-U'it,0) = -g{t) 



f E r(p'(r,0-))dr 



which shows that — [/'(t, 0) can he interpreted as the total quantity of cars passing through the 
junction point X = during the time interval [0,t). As a consequence, the quantity —U^{t,0) can 
also be interpreted as the instantaneous flux of cars passing through the junction point. 

We now give a further interpretation of the problem in the special case m ~ 1. In the special 
case ?7i = 1, imagine for a moment, that we come back to a discrete description of the traffic, 
where each car of label k has a position Xk{t) with the ordering Xk{t) < Xk+i{t) < 0. We can be 
interested in the label k of the car Xk{t) < which is the closest to the junction point X = 0. 
Let us call it K{t). We can normalize the initial data such that K{0) = 0. Then the quantity of 
cars that have passed through the junction point X = during the time interval [0, t) is equal to 
—K{t), which is the exact discrete analogue of the continuous quantity ~U^{t, 0). 

On the other hand the number of cars between the positions a ~ xa(^) and h = xsit) is 
obviously equal to i? — A, and its continuous analogue on the branch ? = ttj = 1 with 7^ — 1, is 
/^ p^(t,X) dX — U^{t,b) — U^{t,a). This shows that U^{t,X) can be interpreted as the exact 
continuous analogue of the discrete labeling of the cars moving in the traffic. 

This interpretation is also meaningful on the "exit" branches, i.e. for j g {m + 1, ...,to + n}. 
Indeed, for such j's, U^{t, X) is the continuous analogue of the discrete label of the car that have 
decided to choose the branch j and which is at time t close to the position X > 0. 

2.2 Getting the Hamilton- Jacobi equations 

We now set 



u^{t,X) = -W{t,-X), X>0, for i = l,...,m 
u^(t,X) ^ ~W{t,X), X>0, for j ==m+l,...,m + n 



and we define the convex Hamiltonians 



Mp) = -^PiYp) 



for i = 1 , . . . , TO 



Hjip) 



ip(_yp) for j ^TO + l,...,TO + n 



(2.; 



(2.9) 



Then we deduce from (2.5) that we have 



Ut+Hkiux)^0, X > 0, for k = l,...,m + n 



(2.10) 



with equality of the functions at the origin, i.e. 



1 



Notice that for the choice Knax 
k G In, that 

ircf(<z) = i(l+9)2 



{t,0) — u{t,0) for any fc e {1, ..., m + n} . 

nax in ^M, we get with f''{p) = f{p) = p{l - p) for aU 



for i = 1, 



L]{q)^ ^LrA-q) for j = m + l,...,TO + n 



In particular this shows that the -/j/c(0) are not all the same, even in the simplest case. 



2.3 The junction condition and its interpretation 



A junction condition is still needed so that the solution of (2.10) be uniquely defined. Indeed, at 
first glance, one may think that Ut{t,0) is equal to —g'{t) which is given by (2.6) (where we have 



assumed (2.7)). The point is that this condition can not be satisfied for every time. One way to 
be convinced oneself of that fact is to consider the case m = n = 1 with /^ = P = /• Then, 
we look at solutions u of the Hamilton-Jacobi equation in M with the artificial junction. We can 
simply associate with it the classical conservation law on the whole real line. We can then consider 
a single shock moving with constant velocity for the conservation law. When this shock will pass 
through the junction point (let us say at time to)j this will mean that Uf (tg, •) is discontinuous in 



space at the junction point. In particular the formal computations of Subsection 2.1 are no longer 
valid at that time tg, even if they are valid for t ^ to. For a general problem, one may expect that 
our formal computations are only valid for almost every time (even if it is not clear for us) . 

In view of Theorem |1.4[ if we restrict our attention to continuous solutions u, then we will 
have u < Uoc where Uoc is the solution associated to the optimal control problem. This shows in 
particular that we have 

U(t,0) <Uocit,0) 

which means (in view of (2.8) and the interpretation of — [/' given in subsection 2.1) that we have 



a universal bound on the total amount of cars passing through the junction point during the time 
interval [0,i). If we assume moreover that this amount of cars is maximal, then we can choose 
(and indeed have to choose) u = Uqc and the natural junction condition is then 

H-{u'),{t,0+))^0 



ut{t,0) 



max 

k = l....,N 



(2.11) 



with 



Hf. (p) = sup {pq ~ Lk{q)) and Lk[p) = sup [pq - Hk{q)) 



q<0 



qei 



Using our assumption (2.1 ) on the functions / , let us define for k — 1, ...,N the Demand functions 



and the Supply functions 



fhip) = 



flip) 



fip) 

fHp'c) 

fHp'c) 
fHp) 



for p < Pc 

for p> Pc 

for p < p^ 

for p> p^- 



From assumption (2.1) on the functions f^, we deduce that 

^fhh'p)^ for i = l,--.,TO 

for j = m+l,...,m + n = N. 



H-{p) 

HJ{p) = -^fs{-l'p), 



Condition (2.11) means that 



[//(t,0)=ut(t,0)= min -H-{u\{t,Q+)) 



fe=i 



N 



= mm I min l/l,(p^(i, 0")), . min J-/^(pJ(t, 0+)) ) . (2.12) 

2— l,...,?n 'y J— ?7i+l,...,m+Ti ^■' 



Notice that from (2.7), it is natural to compare 

1 



J\p\t,Q-)) and p{p>(t,Q-^)). 



Then condition (2.12) is nothing else that the Demand and Supply condition of Lebacque, which 



claims that the passing flux is equal to the minimum between the Demand and the Supply, as it 
is defined in [53] (at least in the case m = 1). 

In the special case m = 1, it is explained in [23] that this condition (2.12) is also equivalent 



to the condition defining the Riemann solver at the junction point in the book of Garavello and 
Piccoli [H]. Let us notice that this condition is also related to the Bardos, Le Roux, Nedelec [S] 
boundary condition. 
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3 Viscosity solutions 



In this section, we consider a weaker assumption than (Al). We introduce the following assump- 
tion: 



(Al') For each i £ In, 

• the function Hi ; M 

• there exists Pq G 
[pI,+oo); 



is continuous and lim 



IpK+c 



H^ip) 



-oo; 



such that Hi is non-increasing on (—00, pp] and non-decreasing on 



When (Al') holds true, the function H- is defined by H^ (p) — infq<o Hi{p + q). We now make 
the following useful remark whose proof is left to the reader. 

Lemma 3.1. Assumption (Al) implies Assum,ption (Al'). 



Next we give equivalent definitions of viscosity solutions for (1.1 1. We give a first definition 



where the junction condition is satisfied in "the classical sense" ; we then prove that it is equivalent 
to impose it in "the generalized sense" . It is essential if one expects solutions to be stable. 

We give a first definition of viscosity solutions for ( |1.1[ ) in terms of test functions by imposing 
the junction condition in the classical sense. We recall the definition of the upper and lower 
semi-continuous envelopes u* and u* of a function u : [0, T) x J: 

u*{t,x)= limsup u{s,y) and u.^,{t,x) = liminf u(s,y). 

Definition 3.2 (Viscosity solutions). A function u: [0, T) x J — > M is a sub-solution (resp. super- 
solution) of (1.1) on Jt if it is upper semi- continuous (resp. lower semi-continuous) and if for 
any (f) G CI{Jt) such that u < (f> in B[P, r) for some P — (i, x) G Jt, r > and such that u = (f> 
at P £ Jt, we have 

4>t{t,x) + H{x,(l>.^{t,x))<{) (resp. > Oj, 



that is to say 

• if X (z J* , then 

• if X — 0, then 



c^t{t,x)+H,{c^^{t,x))<0 



>o;; 



0t(i,O) 



max H^ 



K{t,0)) <0 (resp. >0). 



(3.1) 



A function w: [0,T) x J — > M is a sub-solution (resp. super-solution) of (1.1|-(1.2) on Jt if it 
is a sub-solution (resp. super-solution) of ( |1.1[ ) on Jt and moreover satisfies w(0, •) < uq (resp. 
m(0, •) > uo). 

A function u: [0,r) x J — > M is a (viscosity) solution of (1.1) (resp. ( |1.1[ )-( [L21 )) on Jt if u* is 
a sub-solution and u^, is a super- solution of (1.1) (resp. (l.l)-(1.2)j on Jt. 



As mentioned above, the following proposition is important in order to get discontinuous 
stability results for the viscosity solutions of Definition |3.2[ 



Proposition 3.3 (Equivalence with relaxed junction conditions). Assume (Al'). A function 
u: Jt -^R is a sub-solution (resp. super-solution) of (1.1) on Jt if and only if for any function 



4> G CI{Jt) such that u < (j) in Jt and u = (j) at (i, x) G Jt, 

• if X £ J* , then 

cl>t{t,x) + HdMt,x)) < (resp. > 0) 
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• if X ~ 0, then either there exists i £ In such that 

Mt,0) + H^{(l),{t,0)) < (resp. > 0) 



or (3.1 ) holds true. 



Proof of Proposition \3.3\ We classically reduce to the case where the ball B{P, r) is replaced with 

Jt. 

The "if" part is clear. Let us prove the "only if" one. We distinguish the sub-solution case 
and the super-solution one. We start with super-solutions since it is slightly easier. 

Case 1: super-solution case. We consider a test function e CI{Jt) such that u > (p vn Jt 
and u = (j) &i {Iq^Xq). There is nothing to prove if a;o 7^ so we assume xo = 0. We have to prove 
that (/>t(io,0) + supjgj^ H~ {(f>l.{tQ,0)) > 0. We argue by contradiction and we assume that 

0t(to,O) -I- sup ffr(0;(io,o)) < 0. (3.2) 

Then it is easy to see that there exists a function (j) e CI{Jt) such that 4> > (f) with equality 
at the point (to,0) and such that 

0^(to,O) = min((/);(io,O),pj,) and 0t(to,O) = 0t(io, 0). (3.3) 

Notice that 

H-{^,{t^,Q)) < HSlito^O)) < H-{^litoM- (3.4) 

The first inequality is straightforward. To check the second inequality, we have to distinguish two 
cases. Either we have ^^(io,0) < (/)^(io,0), and then <^^(io,0) = p^ and we use the fact that the 
minimum of H[^ is Hi{p\^). Or (/)^(io,0) = (/!)^(to,0) and then this common value belongs to the 
interval (— (X),Pq] on which we have Hi = Hf. 

Since u > 4> m. Jt and u = (/) at (io, 0), we conclude that either 

0t(io,O)-f supiJ-(^^(io,0))>0 

or there exists i (z In such that 

Mto,0)+H,{4>Uto,0))>0. 



In view of (3.3 1 and (3.4), we obtain a contradiction with (3.2) 



Case 2: sub-solution case. We consider a function <j> G C^{Jt) such that u < (f> in Jt and 

u — (j) ai (to^xo). There is nothing to prove if xq 7^ and we thus assume xq = 0. We have to 
prove that (j)t{to,0) + supj^j^ II~ (0^(io, 0)) < 0. We argue by contradiction and we assume that 

Mto,0) + sup Hr{(l)l(to,0)) > 0. (3.5) 

In order to construct a test function (j), we first consider /i C In the set of j's such that 

iJ7((/)i(to,0))< 8upHr{cbl{to,0)). 

i&In 

Since Hj is coercive, there exists q^ > pg such that IIj{q^) = supjgj^ H~{(j)l.{to, 0)). 

We next consider a test function cp S CI{Jt) such that (p < (j) with equality at {to, 0) and such 
that 
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Notice that for all j E In, 

H,{^i{to,0)) > sup Hr{4>l{to,0)) = sup H-{cj,l{to,0)) (3.7) 

where for the inequality, we have in particular used the fact that Hj is non-decreasing on [pq, +cxd). 
Since u < cj) in Jx and u — (f) at {to,0), we conclude that either 

Mto.O) + snp H-{4>l{to,0)) < 

■iGIn 

or there exists j G In such that 



In view of (3.6) and (3.7), we obtain a contradiction with (3.5). This ends the proof of the 



Proposition. D 

We now prove Theorem |1.4| 



Proof of Theorem \1.4\ Let us consider a function (j) € CI{Jt) such that u < 4> with equality at 
(to,0) with ig e (0,T). Modifying <f) if necessary, we can always assume that the supremum of 
u — (f) is strict (and reached at (ig, 0)). For ry — (771, ..., ry^v) e (M+)^, we set 

Mj, = sup ( u{t, x) — (j)(t, x) — p^ 

{t,x=Xjej)eJT \ Fj 

Because u is continuous at (io,0), we get for 7; e (M+)^ that 

M^ -> A/o = , , , . . n /Q o^ 

, „' „^ , _^ as soon as one 01 the component rji^ — ;> (J. (3.oj 

(^i , a; j y (^tQi t)j 

where (t^,x^) d Jt is a point where the supremum in Af^ is reached. 



Moreover given the components rjj > for j £ In \ {io}i we can use (3.8) in order to find 
r]ig > small enough to ensure that x'' e J* . Then we have in particular the following sub-solution 
viscosity inequality at that point (i'',x''): 

0. + H,o ('^. - ^) < 0. 

Therefore passing to the limit r/i,, — > 0, we get 

c^t + H-{cj,l")<0 at (io,0). 

Because this is true for any ig G In, we finally get the sub-solution viscosity inequality at the 
junction: 

(f>t + uia^Hr{(j>l) < at (tg, 0). 

iEIn 

Now the fact that u < Uqc follows from the comparison principle. This ends the proof of the 
Theorem. D 

4 The minimal action 

We already mentioned that the optimal control solution of the Hamilton- Jacobi equation defined 



by ( 1.8 ) plays a central role in our analysis. We remark that for x e J and t > 0, 



Uoc(i, x) = inf {uoiy) + V{0, y; t, x)} (4.1) 
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where 

rt 



V{0,y;t,x)= min / L(X(r), X(T))rfr. 

XGAiO,y;t,x)jQ 

More generally, keeping in mind the weak KAM theory, we define the so-called minimal action 
V : {{s, y, t, x) e ([0, cx)) x J)^, s < t} ^ M by 

2?(s,2/;i,x)- inf f L{X{T),X{T))dT. (4.2) 

It is convenient to extend 2? to {s = i]. We do so by setting 

V{t,y,t,x) = i 

I +0O II y ^t X. 

Remark 4.1 (Dynamic Programming Principle). Under assumptions (AO)-(Al), it is possible (and 
easy) to prove the following Dynamic Programming Principle: for all a; G J and s € [0,i], 

Uoc{t, x) = inf {uoc(s, v) + 2?(s, y; t, a;)} . 
ye. J 

Notice that a super-optimality principle will be proved in Proposition |5.1| 

Before stating the main result of this section, we 

The following result can be considered as the core of our analysis. The most important part 
of the following theorem is the fact that the minimal action is semi-concave with respect to (t, x) 
(resp. {s,y)). 

Theorem 4.2 (Key inequalities for T)). T) is finite, continuous in {(s, y;t,x): < s < t,x,y d J} 
and lower semi- continuous in {(s, y;t,a;): < s < t,x,y d J}. Moreover, for all (so,2/o) and 
(to,xo) € (0,T) X J, So < to, there exist two functions 0, V' G CI{Jt) and r > such that 

• 4>> 'D{so,yo; •, •) on a ball B{Po,r) with equality at Pq — (tQ,xo) and 

^t+H{xo,^,)>0 at {to,xo). (4.3) 

• ij-) > 'D{-,-]tQ,XQ) on a ball B{QQ,r) with equality at Qq — (so,yo) o.nd 



-A + H{yo,-i;y)<0 at (so^yo) if N>2, 
-^s + Hi{-^y)<0 at {so,yo) if N = I. 



(4.4) 



Moreover, for all R > 0, there exists a constant Cr > such that we have 

d{yo,XQ)<R =^ \(l)x{to,xo)\ + \ipy{so,yQ)\<CR. (4.5) 

Remark 4.3. As we shall see when proving this result, we can even require equalities instead of 



inequalities in (4.3) and (4.4). 



Since the proof of Theorem 4.2 is lengthy and technical, we postpone it until Section [6] When 



proving the main results of our paper in the next section, we also need the following lower bound 
on v. We remark that this bound ensures in particular that it is finite. 

Lemma 4.4. Assume (Al). Then 

2?(s, y; t, x) > ^,^_ g\ d{y. xf - Co{t - s) 

2 

where Cq : — max(0, — Lo(0) + — ), 7 appears in (Al), jq — max^g/^ |i7';(0)| and Lo(0) is chosen 
as in ( l.lOj^ . Moreover, 

V{s,x;t,x) <Lo{0){t~s). 

In particular, if {tn, Xn) -^ (t, x), then 2?(i„, a:„; t, x„) — > as n — > cx). 
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Proof of Lemma \4-4\ We only prove the first inequality since the remaining of the statement is 
elementary. We have 

2 

Hp) > h^ + L'Mp + HO) > Ip" - 70 bl + io(0) > V + io(0) - ^. 



This shows that 



U{p)>\p'-Co. 



Thus we can write for X{-) G A{s, y; t, x), 



L(X(T),X(r))dT>-Co(t-,s) + ^ / (X(T))2dr. 



Then Jensen's inequality allows us to conclude. 



(4.6) 



D 



5 Proofs of the main results 



In this section, we investigate the uniqueness of the solution of (!.!)-( 1.2). In particular, we will 



show that the solution constructed by Perron's method coincide with the function Mqc coming 
from the associated optimal control problem. 

5.1 Super-solutions and super-optimality 

In this subsection, we will show that a super-solution satisfies a super-optimality principle. For 
the sake of clarity, we first give a formal argument to understand this claim. We consider the 
auxiliary function, for s < t. 



Ut.x{s) ^ mi{u{s,y) +'D{s,y;t,x)} 

y€J 



(5.1) 



and we are going to explain formally that it is non-decreasing with respect to s as soon as u is 



a super-solution of (1.1). We call this property a super-optimality principle. Notice that this is 



strongly related to the fact that the quantity Ut^x{s) is constant in s if m is equal to the optimal 
control solution Uqc- 

Assume that the infimum defining U is attained for some y £ J. Then we write 

^t,x(s) = dsu{s,y) + d,V{s,y;t,x) 
dxu{s,y) = -dyV{s,y;t,x). 



Moreover assuming V to be smooth (which is not the case), we formally get from (4.4) the fact 
that ds'D{s,y;t,x) > H{y,—dy'D{s,y;t,x)) (at least in the case N > 2). Hence 

U'is) > dMs, V) + H{y, dxu{s, y)) > 0. 



We thus conclude that Ut^x is non-decreasing if u is a super-solution of (1.1). We now give a 
precise statement and a rigorous proof. 

Proposition 5.1 (Super-optimality of super-solutions). Assume (Al). Let u: [0,r) x J^-M he 

(5.2) 



a super- solution of (1.1) on Jt such that there exists cr > such that for all {t,x) £ Jt, 

u{t,x) > -cr{l + d{x,0)) 



Then for all (t, x) G Jt o,nd s G (0, t], 



u{t,x)> vai{u{s,y) -^T>{s,y;t,xy\ 
ye.] 



(5.3) 



Assume moreover (AO) and that u is a super- solution of (l.l)-[ 1.2) on Jt- Then we have u > Uoc 
on [0,T) X J. 



15 



Proof of Proposition\5. 1\ The proof proceeds in several steps. 



Step 1: preliminary. Notice first that from (6.7), we get 



u{s, y) + V{s, y- 1, x) > J d{y, xf - Coit - s) - a(l + \y\). 

Using the lower semi-continuity of P, we see that the infimum in y of this function is then reached 
for bounded y's. Moreover by lower semi-continuity of the map (s, y; t, x) i— >■ u{s, y) + 'D{s, y; t, x), 
we deduce in particular that the map {s;t,x) i— > Ut,x{s) (and then also s i— ?> Ut^x{s)) is lower 
semi-continuous. 

Step 2: the map s H> Ut,x(s) is non-decreasing. We are going to prove that for s G (0,t), 
Ut xi^) ^ in the viscosity sense. We consider a test function (p touching Ut^x from below at 
s € {0,t). There exists y such that 

Ut,x{s) = u{s, y) + Vis, y; t, x). 

We deduce from the definition of Ut^x that 

(p(s) - V{s, y; t, x) - [ip{s) - V{s, y; t, x)] < u{s, y) - u{s, y). 



By Theorem 4.2 there exists a test function -0 such that il) > !?(•, -^tjx) on a ball B{Q,r) with 
equality at Q = (s, y). Hence, we can rewrite the previous inequality by replacing V with tfj. We 
then obtain that (s^y) i— >■ (p{s) — il){s,y) is a test function touching u at (s, y) from below. Since 



M is a super-solution of (1.1 ), we have in the cases A^ > 2 or A^ = 1 and y 7^ 

/(s) > ^s{s, y) - H{y, ~dy^{s, y)) > 
and in the case N = 1 and y = 

^'{s) > A{s,y) - H^{~dyi:{s,y)) > Vs(s,y) - H^{-dy^s,y)) > 
where we used the properties of the function i/j given by Theorem |4.2[ 



Step 3: conclusion. Let us define for (i, x) e Jt the following kind of lower semi-continuous 
envelope (for the past in time) 

u{t,x) = liminf{u(t„,a;„) : (i„,a;„) -> {t,x),tn < t}. 

Let us notice that we have 

u^, — u^, — u on Jt- (5-4) 

Given a point (i,x) € Jt, let us consider a sequence (i„,a;„) -^ {t,x) such that 

u(t,x) = lim u{tn,Xn)- 

n— 5-+00 

Using Lemma |4.4[ we have for any s < t^ < t 

Ut,x^{s) < Ut^xA^n) < u{tn,Xn) + X>(t„ , a;„, i, x„) -^u{t,x). 

Therefore from the lower semi-continuity of U, we get 

Ut,x{s) < u{t,x). 

Again from the lower semi-continuity of the map {t,x) 1— > Ut^x{s), we get passing to the lower 
semi-continuous envelopes in {t,x): 

Ut,x{s) < u^{t,x) ^ u{t,x) 



where we have used (5.4). This shows (5.3) for < s < t. This is still true for s = t by definition 



of v. The proof is now complete. D 
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5.2 Comparison with sub solutions 



Proposition 5.2 (Comparison with sub solutions). Letu: Jt ^- ^ be a sub-solution of (l.l\-{1.2) 
on Jt, such that there exists cr > such that for all {t, x) G Jt, 



u{t,x) <cr(l + d(a;,0)). 
Then we have u < Uqc on Jt- 

In order to prove Proposition |5.2[ we first state and prove two lemmas. 



(5.5) 



Lemma 5.3. Assume (AO)-(Al). Then the function Uqc defined in (1.8) satisfies 

\uoc{t-,x) - uo{x)\ < Ct. 



Proof of Leninia \5.3\ We first get a bound from below. Using (6.7), we deduce (denoting by Lug 
the Lipschitz constant for uq): 

uo{y) +'^{s,y;t,x) > uq{x) + ^Ad{y,x)f -Cot- Luod{y,x) 



with 



This implies that 



-c. 



inf 

aG[0,+oo) 



> Uo{x) - C2t 

\-a^ - Co - Luaa\ 



> — oo. 



Uoc{x) > Uo{x) - C2t. 

We next get a bound from above. We have 

Woc(a^) < "0(2;) + T){0, x; t, x) < uo{x) + Mt 

with 

M ^ sup L,(0). 

ieiN 

This ends the proof of the lemma. 

Lemma 5.4. Assume (AO)-(Al). Let u : [O.T) x J - 



D 



be a sub-solution of (1.1|-(_Z.^) on J' 



satisfying (5.5 I. Then there exists a constant C > such that 

u{t,x) < uq{x) + Ct for all {t,x) E Jt- (5-6) 

Proof of Lemma \5.4\ Using the Lipschitz regularity of uq, we can easily consider a smooth ap- 
proximation Uq of Uq such that Mq > Uq and |mq — uo|ioc(j) — > as 77 -> 0. Then consider the 
following supremum for r],a > 



Nri,a= sup {u{t,x) — u^{x) — Ct — ad{x,QY 

{t,x)e[Q,T)xJ 



V 



T-t 



}• 



We claim that Nj^^ ^ for some C large enough independent on ry, a > small enough. The 
lemma will be obtained by letting a and ry go to 0. We argue by contradiction and assume that 



Nn,a > 0. Thanks to (5.51, the supremum iVj,^^ is attained for some [t^x) S [0,T) x J. If t = 0, 



we have N^^a < 0. Therefore t > and we can use the fact that u is a sub-solution to obtain for 

X — Xj^Ci 

^ + C - umxL,{0) < ^ + C + H{x,d^ul{x) + 2ax,) < 

where we have used assumption (Al) to estimate H from below. Notice that we have also made 
use of a slight abuse of notation in the case a; = 0. Choosing C = maxjg/^ |ij(0)| allows us to 
conclude to a contradiction. This ends the proof of Lemma 5.4 D 
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We now turn to the proof of Proposition |5.2[ 

Proof of Proposition \5.S\ The proof proceeds in several steps. 
Step 1: preliminaries. Let us consider 

M — sup {u(i, x) — Uoc{t, x)} . 

(t,x)e[0,T)xJ 

From Lemmas |5.3| and |5.4[ we deduce that we have M < 2CT < +00. We want to prove that 
M< 0. 

To this end, we perforin the usual corrections considering the following supremum for ?/, a > 



Mri,a ~ sup <.u{t,x) — Uoc{t,x) ~ ad{x,0) — 



{t,x)e[o,T)xj L T - t 

As it is proved classically, we also have that A/^^q, — ^ -M^.o as a — > where 

Mrffi = sup <^ u(i, x) - Uoc{t, x) 



{t,x)elo,T)xj L T - t 

We argue by contradiction by assuming that Af > and then Af^^o > M/2 > for 77 > small 
enough and fixed for the rest of the proof. 

Step 2: reduction to t > 0. Notice that the supremum M^j^a is achieved for points (t, x) € 



[0,r) X J. Using again Lemmas 5.3 and 5.4 we also deduce that 



M/2 < A/^,0 < Mrj^a + Oa{l) < 2Ct 

and hence t > ^ > for a small enough. 

Step 3: a priori bounds. Using the argument of Step 1 of the proof of Proposition [O] we see 
that there exists y G J such that 

Uoc{i,x) = uo{y) +V{0,y;i,x). 

Therefore we can rewrite Afjj.a as 

71 
Mjj,a= sup {u{t,x)-uo{y)-'D{0,y;t,x)-ad{x,0)'^ - - — -}. 

0<t<T,x,yeJ 1 — t 

and the supremum is achieved for (t, x, y) G (0, T) x J^. Notice that this supremum looks like the 
classical one for proving the comparison principle for viscosity solutions, with the usual penaliza- 
tion term (y — x)'^ /e replaced here by the function X'(0, y; t, x). 



In view of the bound (6.7) from below on T) and (5.6), we derive from Af^ q, > that 



"^ + ad{x, Qf + Jd(y, xf < Cot + Ct + Lu„d{y, x) 



T-t ' ' ' At 

where L^q denotes the Lipschitz constant of uq. We conclude that there exists Ct such that 

ad{x,Of<CT and d{y,x)<CT (5.7) 

where Ct depends on T, Co,C, L„„ and 7. 
Step 4: getting the viscosity inequality. Since t > 0, we have in particular that 



u{t,x)- [V{Q,y]t,x)+ad{x,Of + ^ 



T-t 



< uil x) ~ [ I?(0, y; i, x) + ad{x, 0)^ + ^^ 
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By Theorem 4.2 there exists a test function (f) such that cj) > I'(0,27; •, •) on a ball B{P,r) with 
equality at P = {i,x). Hence, we can rewrite the previous inequality by replacing V with (j). We 
then obtain that {t, x) i-)- (f){t, x) + ad{x, 0)^ + ^^ touches u from above at {t, x) with f > 0. We 
use next that w is a sub-solution of ( 1.1 ) and get for x = XiCi 

V 



y2 



+ Mt, x) + H{x, (j)^{t, x) + 2ax.,) < 



where we have made use of a slight abuse of notation in the case a; = 0. On the other hand, we 
have 

Mt,x) + H{x,(l>^{i,x))>0 

therefore 



y2 



+ H{x, (t>^{t, x) + 2axi) - H{x, (t>^{t, x)) < 0. 



On the one hand, from (5.7), we have < axi < ^faCr- On the other hand, we can use (5.7) and 



(4.5) in order to conclude that 



|(/.,(i,x)|<C 

for some constant C which does not depend on a. We can now use the fact that the Hamiltonians 
are locally Lipschitz continuous in order to get the desired contradiction for a small enough. This 
ends the proof of the Proposition. D 



5.3 Proof of the main results 

In this subsection, we prove the main results announced in the introduction. 



Proof o] Theorem 1.2. We simply apply Propositions 5.1 and 5.2 and get u < Uqc < v which 
implies the result. D 

In order to prove Theorem |1.1[ we should first prove that solutions are Lipschitz continuous. 



Lemma 5.5. Assume (AO)-(Al). Let u be a solution of (1.1)-(1.2) on Jt- Then u is Lipschitz 
continuous with respect to (t,x) on Jt- 



Proof of Lemma\5.5\ We first recall (see Lemma 3.1) that (Al) implies (Al'). We know that the 



solution u = u* given by Theorem A. 2 satisfies for some constant C > and all (t,x) g Jt, 



\u{t,x)-uo{x)\ < Ct. 

From the comparison principle (Theorem |1.2[ ), we deduce that u — u* < u^, and then the solution 
u — u* ^ u^, is continuous. 

For Hq > small (with Hq < T), we now consider h £ (0, Hq) and 

v{t, x) = u{t + h,x) — sup(u(ft-, x) — uo{x)). 



This new function satisfies in particular v{0, x) < uq{x). Therefore u is a sub-solution of ( 1.1 )-( 1.2 1 
on Jx-ha- We thus conclude from the comparison principle that w(i,x) < u{t,x), which implies 

u[t + h,x)<u{t,x) + Ch for all {t,x) £ Jr-ho- 

Arguing similarly, we can prove that u{t+h, x) > u{t, x) — Ch. Because Hq can be chosen arbitrarily 
small, we conclude that u is Lipschitz continuous with respect to time on the whole Jt- 



Since m is a viscosity solution of (1.1), it satisfies in particular (in the viscosity sense) for each 
i e In- 

H,{u^)<C on (0,r)xJ*. 

This implies that there exists a constant C such that (in the viscosity sense) 

\ux\<C on (0,r)xJ*. 

This implies that u is Lipschitz continuous with respect to the space variable. This ends the proof 
of the lemma. D 
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We now turn to the proof of Theorem |1.1[ 

Proof of Theorem \l . 1\ The uniqueness of the solution follows from Theorem |1.2| The existence is 
obtained thanks to the optimal control interpretation (uqc is a solution). The Lipschitz regularity 
was proved in Lemma 5.5 above. The proof of Theorem |1.1| is now complete. D 



Proof of Theorem 1.3. The fact that the solution is equal to Uqc follows from Propositions 5.1 and 



5.2 The representation formula (1.9 1 follows from (4.1). D 



6 A complete study of the miniinal action 

6.1 Reduction of the study 

We start this section with the following remark: the analysis can be reduced to the case {s,t) = 
(0, 1). Precisely, using the fact that the Hamiltonian does not depend on time and is positively 
homogeneous with respect to the state, the reader can check that a change of variables in time 
yields the following 

Lemma 6.1. For all x.,y € J and s <t, we have 

V{s, y; t, x) = {t- s)V U, -^; 1, -^"j . (6.1) 

\ t — s t — s J 

This is the reason why we consider the reduced minimal action 2?o : J'^ — >• K defined by 

Vo{y,x)^V{Q,y-l,x). 

Thanks to the previous observation, it is enough to prove the following theorem in order to get 
Theorem 14.21 



Theorem 6.2 (Key equalities for "Dq). Let us assume (Al). Then for all y,x Cz J, the 'Dfj{y,x) 
is finite. It is continuous in J^ and for all y,x C^ J , there exists a function ipQ G Cl{J'^) such that 
</'o ^ 2?o on J- , ipoi^y, x) — I'o(j/, 2;) and we have 

ifx^O: {ipo - xdxipo - ydyipa){y,x) + H{x,da;(po{y,x)) = (6.2) 

ifx^O: {fo- xdxipo-ydyipn){y,0)+ sup H~{dl.ipQ{y,0)) ^0 (6.3) 

ieiN 

and if y ^ 0, 

{ifo - xda:(po - ydy(pa){y,x) + H{y,~dyLpo{y,x)) = (6.4) 

and if y = 0, 



{ipo - xda;ipQ - ydy<fio){Q,x) + sup^^j^ H^ {-dlipo{0,x)) =0 if N >2, 
(ipo - xdxfo - ydy(po){0, x) + ^i(-9^¥'o(0, a;)) =0 if N = 1. 



(6.5) 



Moreover, for all R > 0, there exists C'n > such that for all x,y €1 J , 

d{y,x) < R^ \dxipo{y,x)\ + \dyipo{y,x)\ < Cr. (6.6) 

Remark 6.3. If Jo — In, then we have Vq G Cl{J'^). This good case corresponds to the case where 
all the Li(0)'s are equal. 

We can interpret Lemma [44] as follows. 
Lemma 6.4. Assume (Al). Then 

Vo{y\x)>'^-d{y,xf-Co (6.7) 

where constants are made precise in Lemma \4.4\ 
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6.2 Piecewise linear trajectories 

We are going to see that the infimum defining the minimal action can be computed among piecewise 
Hnear trajectories, and more precisely among trajectories that are linear as long as they do not 
reach the junction point. This is a consequence of the fact that the Hamiltonians do not depend 
on X and are convex (through Jensen's inequality). 

In order to state a precise statement, we first introduce that optimal curves are of two types: 
either they reach the junction point, or they stay in a branch and are straight lines. This is the 
reason why we introduce first the action associated with straight line trajectories 

U{x,-y,) if (2/,x)e J2\{(0,0)}, 
I'straight(2/,a;) = ■^ Lo(0) if y = Q = x, 

^oo otherwise 

and the action associated with piecewise linear trajectories passing through the junction point 



2?junction(y,a;) = inf {fi(Ti,2/) + '?2(t2, a;)} 

0<Ti<r2<l 



where 




and 



f^ j - riLo(O) for y = y,e, ^ 0, n ^ 
for y = 
for Ti = 0, y 7^ 

(1 - T2)L, (j^) +T2LoiO), for X = x,e, ^0,T2^1 
f2(T2,a;) = <^ i„(0) forx^O 

for T2 — l,x =/: 0. 

Remark 6.5. By defining the f^'s in such a way, we treat the degenerate cases: x = or y = 0. 
Indeed, ti (resp. T2) measures how long it takes to the trajectory to reach the junction point 
(resp. the final point x) from the starting point y (resp. the junction point). 

The following facts will be used several times. 
Lemma 6.6. The function £1 (resp. £2) is continuous in (0, 1] x J* (resp. in [0, 1) x J* ). 
Lemma 6.7. The function £i, i — 1,2 are lower semi- continuous in [0, 1] x J . 
Proof. Consider the function defined for (t, y) G [0, 1] x J by 

r rL,(-^) ify = y,:e, ^0,r^0 
g{T,y)^l tLo{0) if y = 

[ +CX) if y 7^0,T = 0. 

From the inequality for r > (consequence of ( |4.6[ )): 

g{T,y) >- Cot, 

4 T 

we deduce that g is lower semi-continuous. Consequently, the map £1 is lower semi-continuous. 
We proceed similarly for £2- D 

We first show the main lemma of this subsection. 

Lemma 6.8. The infimum defining the reduced minimal action Vq can be computed among piece- 
wise linear trajectories; more precisely, for all x,y € .J, 

X'o(y, x) = min(I?straight(y, a;), Injunction (y, a;)) ■ 
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Proof. Wc write with obvious notation Vod/, x) = 'i'n-ixeAoiy,x) ^i^)- I^^ order to prove the lemma, 
it is enough to consider a curve X e AQ{y,x) and prove that 

£{X) > min(I?straight(y,a;),I?junction(y,a;))- 
To do so, we first remark that the uniform convexity of Li imphes that for all po G I^i we have 

U{p) > L,{po) + L',{po){p -po) + l{p- Pof- (6.8) 

We now consider an admissible trajectory X: [0, 1] — )• J and we treat different cases. 

Case A: X((ti,i2)) C J*. We assume that a curve X stays in one of the branch J* on the time 
interval (ti,i2) with ti < ^2- In such a case, we consider the curve X with same end points X{ti) 
and -'^'(^2) in Ji but linear. If pg G M is such that PqCi — X{t) for t G (^1,^2) and pci — X{t), we 



deduce from (6.8 1 that 



' L{X{T),X{T))dT > [ ' L{X{t), X{T))dT +1 r \X{t) ~ X{t)\Ht. (6.9) 

ti Jti ^ Jti 

Case B: X{[ti,t2]) C Ji with X(ti) = X(t2) = 0. In that case, let usset X(t) = for r g [^1,^2]- 
Using (6.8 ) with po — and the definition of Lq as a minimum of the Lj's (see ( 1.10 )), we get that 



Hp)>Lo{0) + mO)p+^P^ 



from what we deduce that (6.91 still holds true. 



Case C: the general case. By assumption, we have X e AQ{y;x) C C([0, 1]). We then 



distinguish two cases. Either ^ -^([0, 1]), and then we define X as in Case A. In this case, (6.9| 
implies that 

£{X) > 'Dstrmght{y,x). 

Or € X([0, 1]), and then we call [ri,r2] C [s,t] the largest interval such that X{ti) = = X{t2), 
and define X as follows: it is linear between and ti, and reaches at ri; it stays at in (ti, T2); 
then it is linear in (r2, 1) and reaches a; at i = 1. Using again the continuity of X, we can find a 
decomposition of [ri, T2] as a disjoint union of intervals Zfe (with an at most countable union) 

[t1,T2] = IJifc 
k 

such that for each k, X{Xk) C Ji^ for some ik € In and X = on dXk- Using Case A or Case B 
on each segment I^, we deduce that 

£{X) > Pj unction (y, a;). 

n 

6.3 Continuity of the (reduced) minimal action 

Lemma 6.9 (Continuity of Injunction)- The function Injunction ** continuous in J^ . 



Proof. We first prove that ^junction is lower semi-continuous. We know from Lemma 6.7 that the 
function 

G{Ti,T2;y;x) ^£i{Ti,y) +£2(t2,x) 

is lower semi-continuous for y,x Cz J and < ^i < T2 < 1. Therefore the function 

Injunction {y;x) = iuf G'(ti , T2 ; J/ , x) 

0<ri<T2<l 
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is also lower semi-continuous (since the infimum is taken over a compact set). Besides, the infimum 
is in fact a minimum. 

We now prove that Injunction is upper semi-continuous at any point {y,x). Consider first 
(ti, T2) e [0, 1]^ such that 

Injunction (y, a;) =fi(ri,y) +£2iT2,x). 

Given any sequence {y^,x^) — >■ {y,x), we want to show that 



We use 



with an appropriate choice of (t{°, Tj ) 



I^junctionly ,X ) < X'junction(y, 2;) + Ofc (1) 



Injunction (/,a;'=) < £i(rf,/) +£2(r2",x'=) 

k ^k\ 



(6.10) 



Case 1: y e J*, x G J*. In this case, we choose {ti,T2) — {ti,T2) G (0,1) and we use 
Lemma |6.6| in order to get 

£i(r^/)^fi(ri,y) 



and 



£2(r2^x")^f2(r2,.T) 



Hence we conclude that (6.10) holds true. 



Case 2: y = 0, X e J* . We choose (rf , t|) = (|y''|, max(r2, \y''\)) € [0, 1)^. We still have t| -^ T2 
and we can use Lemma |6.6| in order to get 



£2iTlx')^£2iT2,x) 



We also have (if y'^ e Jj ) 



£l(r^/)<|y'=|L, 



|/|Lo(0)^0 = fi(ri,0). 



(6.11) 



Hence we conclude that (6.10) holds true. 



Case 3: 2/ G J*, x ^ 0. We choose {t^,t^) = (min(Ti, 1 - |x'=|), 1 - [x'^D G (0, 1]^ We stiU have 
Ti — > Ti and then 

fl(r^/)^fl(rl,y) 

(since f i is continuous in (0, 1] x J*). We also have (if x'^ G Ji) 



f2(T-2^x'=)<|x'=|L, 



+ (1 - |x'=|)Lo(0) ^ Lo(0) - £2(r2,0) 



(6.12) 



Hence we conclude that (6.10) holds true. 



Case 4: y = 0, x ^ 0. We choose (t^.t^) = i\y''\, 1 - \x''\) G [0, 1) x (0, 1]. We deduce ( [6l0 l 
from (16. Ill) and (|6.12|). D 



Lemma 6.10. The function T>q is continuous in J^ 



Proof. Since mstraigiit is lower semi-continuous, we can use Lemmas 6.8 and 6.9 in order to conclude 
that Vq is lower semi-continuous. 

Consider {y, x) G d{Ji x Ji) \ {(0, 0)}. Then either a: = or y = 0. Moreover for y — yiCi and 

V- .■ f,, rW / ^i^^'^)+^2(l,a;) if-T^«=0 \ < r /^._„,^ 

injunction 12/, Xj S <y £^ (Q, y) + f 3 (Q, x) if y, = J " ^ ^ ' ^'^ ' 
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Therefore for each i e /jv, we have for (2/,a;) G d{Ji x J^), 

injunction (y, a;) < 2?straight(y,a;) 
Therefore we have with y = yiei, x — XiCi 



'Dq{jj,x) = < 



min(I?junction(2;, a;),Li(xi - yi)) if {y,x) e JiX Ji 

Injunction (y, a;) if (.y,x) £ d{Ji X Ji) (6.13) 

Injunction (y, a;) otherwise. 



This imphes that I?o is continuous in J^. D 

6.4 Study of Injunction 

In view of ( |6.13[ ), we see that the study of I?o can now be reduced to the study of Injunction- The 
function ^junction is defined as a minimum over a triangle {{ti,T2) € [0,1]^ : ti < T2}. We will 
see below that ^junction is defined implicitly when the constraint ri < T2 is active (I?impiicit) or 
defined explicitly if not (I?iincar)- In other words, it will be linear "as long as" trajectories stay 
some time {t2 — ti > 0) at the junction point. 
We first define for (y, a;) G J^, 

implicit (2/, a:) = inf {£i{T,y) +£2{t,x)} . (6.14) 

0<r<l 

The continuity of I^impiicit will be used later on. 

Lemma 6.11 (Continuity of I?;'' jj^j^). The restrictions I^f^pii^it 0/ ^implicit o-re continuous in 
(Jj X Ji) \ {(0,0)} and continuous at (0,0) if j G /o or i £ Iq. 



Proof of Lemma \6.11\ From Lemma 



6.7 



we deduce that 'C'l^^^ii^n^ is lower semi-continuous on 



Jj X Jj. We now show that T'j^ ^^^-^^ is upper semi-continuous at any point (y, x) e (J, x Ji)\{(0, 0)} 
and also at (0, 0) if j £ Iq or i E Iq. We first consider r e [0, 1] such that 

^implicit(2/.^)=^l(^'2/)+^2(^'^) ^ith 0<T<1. 

For any sequence {y'',x'^) — >■ {y,x) with [y'^.x'^) G Jj x Jj, we want to show that 

^S>piicit(/,^') <^i^;piicit(y,^) + Ofc(l) (6.15) 

Arguing as in Lemma |6.9[ we use 

»fmplicit(y^^')<^l(^^y')+^2(r^x'=) 

and we choose t^ as follows 

iiy e J*,x e J* : r''=TG(0, 1), 

iiy = 0,xeJ:: r'= = |yfe|G[0,l), 

ifyGj;,a; = 0: t*^ = 1 - I-t*^] G (0, 1], 

ifa: = 0,iG/o: t*^ = 1 - |.tH G (0, 1], 

iiy = 0,x = 0,ieIo: t'= = |y'=| G [0, 1). 

This ends the proof of the lemma. D 

We next define for (y, x) G J^ x J* 

^Lar(y>^) = -m^i)y + L'Mt)x + io(0) (6.16) 
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where ^; are defined thanks to the following function (for I £ In) 



Precisely, f; = {^1) ^(0) 7^ when I ^ Iq (sec Lemma 6.17 below). We will sec that Ki plays an 



important role in the analysis of Injunction- In particular, it allows us to define, when i ^ Iq and 
j ^ If), the following convex subset (triangle) of Jj x J,: 

It is convenient to set A^' — 9 ii i E Iq 01 j E Iq- We next state a series of lemmas before proving 
them. 

Lemma 6.12 (Link between ^junction, Ainoar, Ampiicit)- 

^Lction(y, ^) - 1 3-^;-^^. i )!' :^ : n ., . ^ ^ ^ ,. (6-i7) 



'junctionvy,-^; 1 I^fiplicit fe, ^) */ (y, x) G ( J, X J,) \ A^\ 



implicit 



Lemma 6.13 (The equations in the interior). The functions 'D'll^ction' -^Uncar '^^'^ -^i 
convex and C^ in J* x J* anrf, i/2? is one of them, it satisfies for {y,x) E J* x J* 

Viy, x) - xd^Viy, x) - ydyViy, x) + H,{dJ){y, x)) = 0, ,g ^^. 

■D{y,x)^xd,T){y,x)~ydyT){y,x) + H,{~dyT){y,x))=0. ^' ' 

Lemma 6.14 (Study of Pimpiicit)- For {y,x) E J* x J*, there exists a unique t = T{y,x) E (0, 1) 
such that 

'^Lpiicitiy^x) ^^iiT,y)+£2{T,x); 
Moreover, 

f 5.^fmpiicit(2/>^) - L', (e.) w^th e. = T3if^, 



i '9!/^fmplicit(y'^)==-^j(^1') ™^^ ^V^-Tit^}- 

Lemma 6.15 (Study of T). For {y, x) E {Jj x Ji) \ {(0, 0)}, there is a unique t = T{y, x) E [0, 1] 
such that 

Moreover T E C{Jj x J, \ {(0, 0)}) and 

'max(0,l-|^) tf iy,x)Ei{0}xJ*)\Ai 
T{y,x) = ' 



min(l,-J^j tf {y,x)EiJ;x{0})\Ai\ 
Lemma 6.16 (I^function ^^ ^^^ boundary). Then we have I'j'unt.tion ^ ^^i^j ^ '^i) ™^^ 



"9a; ^junction (2/' ^) — ^'i i^^) /g ^g) 

9y'^iunctioniy^x)^-m^y) 



where S.y < < ^x satisfy 



•(max(a;,e+), (if")"! (if. (?.))) J/(y,x) € ({0} x J,) \ A^'^ 

(?.,C.) = <!^^-fl ./(y,x)e({0}xJ.)nA- 

V?.,?.^ ^((if+)-i(X^.(e.)),-max(y,-e7)) */(y,x)G(J, x{0})\A^' 

.(^+,^7) z/(y,x)e(J, x{0})nAJ-\ 
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Moreover we have 

f ^^/LtionlO, x)^f- (L,(^,) - Lo(0)) + Lo(0) for x e J* 

1 2?]'Ltio„(2/:0) = -^(L,(e,)-io(0)) + io(0) for y e J* 

and 

^/unction (^' y) - ^^^r^/unctionCa;^ V) " J^^y^j^nctionCa^- y) 

Lo(0)+if,(max(a;,e+)) */ (y, x) G {0} x J„ 

Lo(0)+if,(-max(y,-e7)) i/ {y,x) e J, x {0} . 



(6.21) 



(6.22) 



Before proving these lemmas, the reader can check the following useful properties of the function 
Ki that will be used in their proofs. 

Lemma 6.17 (Properties of Ki). Assume (Al). Then for any I € In, we have 

K'i{0>l\^\ for ee(-oo,0), 
Ki{0<~l\^\ for ee(0,+(X)). 

We define (Kf)^^ as the inverse of the function Ki restricted to (— oo,0], and {K^)^^ as the 
inverse of the function Ki restricted to [0, +oo). We set 

Then we have 

±ef = if lelo, 
±ef > if lelN\Io. 

Moreover we have 

KiiO = -Hi{L[iO)-LoiO)- (6.23) 

Proof of Lemmas[KJWKJ^ The proof proceeds in several steps. 
Step 1: first study of 'D--^^^^^^^^. Let us define 

G{Ti,T2,y,x) =£i{Ti,y) +£2{t2,x). 

For Ti,T2 G (0, 1), and setting 

r 1 1 - T2 

and Vy = {^y, 0, 1, 0) and Vx = (0, S,x, 0, 1), we compute 



6. = -!^, e. = r^ (6-24) 



D'Gir„r2,y,x) ^ ^^^^V^Vy + ^^VjVx > 0. 

Ti " 1 - T2 

Therefore G is in particular convex on (0, 1) x (0, 1) x J* x J* . Because G is in particular lower 
semi-continuous on [0, 1] x [0, 1] x J* x J*, and 

G{0,T2,y,x) = +oo = G{Ti,l,y,x) for {y,x) e J* x J* , (6.25) 

we deduce that 

^j'unction(y' ^) = „^ inf ^, G(ti, ra, y, x) for (y, x) G J* x J*. 

■' 0<ri<T2<l 

This implies that 'Di^^ction i^ ^^^^ convex in J* x J* . Notice that in particular 

Dl^,^Gin,T2,y,x) = y^L'^iCy) > 
~i 
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and 

The map (ri, T2) i— > G'(ti, r2, y, x) is then strictly convex on the convex set 

{(Tl,T2)e(0,l)^ T,<T2}. 



Therefore using again (6.25) and the lower semi-continuity of G, we deduce that it has a unique 



minimum that we denote by (ti,T2) satisfying < ti < T2 < 1. 
Step 2: study of T>f^ i:^^^^. Let us consider the following function 

e{T,y,x) = G(T,T,y,x). 

For T € (0, 1), setting 

y X 



7 ^jj -1 

r 1 — T 

and proceeding similarly as in Step 1, we can deduce that 

^fmpiicit(2/. 2;) = ^Jnf^^ e(T, y, x) for {y, x) e J* x J* 

and that 2?i'j^jpii(.it is also convex on J* x J*. We can also deduce that the map r i-)- e{T,y,x) is 
strictly convex on (0, 1) for (y, x) e J* x J* and that it has a unique minimum that we denote by 
T £ (0, 1) such that 

Using the derivative with respect to r, we see that r is characterized by the equation 

F = with FiT,y,x):=K,{~^)-K,i-^). (6.26) 

r i — T 

Moreover 

9^i^(r, y, x) = i:'?^e(r, y, x) > 0. 

Using the regularity (7^ of L; given in assumption (Al), we see that the unique solution r — T{y, x) 
of F(T,y,x) = is continuously difFerentiable with respect to {y,x). Therefore we deduce that 
2?^pnciteCi(j;xJ-). 



We have 



^implicit (y'^) = £i{T{y,x),y)+£2{T{y,x),x), (6.27) 

5y2?f;p,icit(y>^) - {dy£^){T{y,x),y) ^ -L'^i^y), (6.28) 

9:^^fmpiicit(y>^) = (9.f2)(r(y,x),x)=iUC.)- (6.29) 



Writing r for T{y^x), and using the optimality condition (6.26), we get 



(^implicit ^^aji'inipiicit Z/'^y^implicit) (2^' ^) 

= rif,(-^) + (1 - r)X,(-^) + io(0) 

T 1 — T 

= if,(-^) + Lo(o) = -i/,(i;.(-^)) 

^ -H,{-dyVt^,,^^{y,x)) 



md.vz.ucjy^^)) 
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j; X J* 



where we have used (6.23) m the second and m the fourth hue. Hence I'ijnpiicit satisfies (6.18) on 



V' 



junction" 



Step 3: further study of V^ 

we have 

^lunction (2/: 2:) = £l (ti , y) + f 2 (t2 , x) 

with < Ti < T2 < 1. Then we can distinguish two cases. 
Case 1: ri < T2. In that case this imphes that 



We concluded at the end of Step 1 that for (?/, x) G J* x J* 



which can be written as 



dri£l{Ti,y) =0, 9r2^2(T2,a;) 







(6.30) 



with ^y < < Cx defin ed in ( |6.24[ ). 

Using Lemma 6.17[ we conclude that (6.30) holds true if and only if Kj{0) > and Ki{0) > 0; 
i.e. j, i £ In\Io- In this case we have ^y = ^~ and S^^ = £,i^ and then 



n 



C- 



T2 = l- 



c 

Sz 



(6.31) 



Moreover, we have in this case Vj^n^^■^n{y, x) = Vf.'^^^^iy, x). 

Using Legendre-Fenchel's equality together with Kj{^j) = and Ki{^^) = 0, we have 



^Unoar(2/' ^) " Vdy^L^AV , x) - xd:,Vl^^^^^{y , 



io(0), 



(6.32) 



and 






Hence V'^^^^^ satisfies ( |6.18[ ) on J* x J* 

Finally we deduce 
A^* n ( J*)^ and then by continuity of I'j'unction ^'^^ ^Uncar' '^^ S'^t 



rom (6.31) that the condition: < ri < T2 < 1 is equivalent to {y,x) £ 



V'"- ,. ^Vi. on AJ\ 

junction linear 



Case 2: ri = T2. If for {y,x) e J* x J* we have 



junction 



(y, x) = f 1 (ti , y) + £2 {t2 , x) 



with Ti = r2, then we have seen that {y,x) £ {J* x J*) \ A^' and Vj^nction\y ^ -^ J - -^implicit 



From Lemma 



6.11 



we also have that 2?fi* li^it ^ ^{Jj ^ -^i) if j £ lo or i £ Iq and in that case 
On the other hand, we have X'/^pii^it ^ C{{Jj x Ji) \ {(0,0)}) if j,i £ In\Io with 
{(0, 0)} £ A'-' in that case. Therefore in all cases we have 

2?^ipUciteC((J,xJ,)\A^'). 

Now from the continuity of Injunction, we deduce that 

^j^nction=2?^iplicit on {J, X J,) \ A^\ 

Step 4: on the boundary (dA^^) n (J*)^- We already know that Injunction is continuous, 
therefore if j, i £ In \Io'- 



linear 



^fmpiicit on 



(y,x) £ Jj X Ji. 



J^ y_ 
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On the other hand, recah that (y, x) E J* x J* , the real r e (0, 1) is characterized by (6.26), i.e 



K,{^y)=K,i^,) with ^y = -^, e. = 



1 



(6.33) 



Notice that if we choose 



we deduce from ^ — 7^ = 1 that 



^y = (j and (^ = Ci' 



(6.34) 



which are o bvious solutions of (6.33). The refor e we conclude that this is the solution. Using 



(|6.28[)-(|6.29[) and the expression of V^^^^^^, (6.34) implies the equality of the gradients of Vf^^^^^. 
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and X'fj'^piieit on the boundary (^A^') n {J*f. Finally this shows that X'j^nction ^ ^'^(J* x J*). 
This ends the proof of the lemmas. 

Proof of Lemma \6.15\ The proof proceeds in several steps. 
Continuity of T. We set for (t, y, x) e [0, 1] x Jj x J^ 

e{T,y,x) =£i(r,j/) +£2{t,x). 



From Proposition 6.12 we already know that there exists a unique t e [0, 1] such that 

^frnpiicitC^' ^) = e(T, y, x) if (2/, x) e J* X J*. 
On the other hand, we have 

f (1 - r)L, (^) + rLo(O) if (y, x) G {0} x J* (case 1), 
\ ri, (-^) + (1 - r)Lo(O) if (y, a:) G j; x {0} (case 2). 

Notice that in Cases 1 and 2, there is a unique r G [0,1] such that 

^i'mpiicit(2/'2;) = e(T,y,x) 



(6.35) 



(6.36) 



and T e [0, 1) in case 1, r g (0, 1] in case 2. Then the continuity of r = T{y,x) in (Jj x Ji)\{(0,0)} 
follows from the lower semi-continuity of e on [0, 1] x Jj x Ji and the uniqueness of t such that 
(lelel) holds. 



Computation of T. We distinguish cases. 

Case 1: {y,x) £ ({0} x J*) \ A^\ Notice that we have 

dre{T,0,x) ^ -Ki{^x) with ^^ 



1 



If a; > C+, then are(T, 0, a;) > and T(0, x) = 0. 

If x < ^,+ , then ^^ = 4+ is a solution of 5re(r, 0, x) = -K»(C:,) = and r(0, x)^l- -^. 

Case 2: (y,x) G (J* x {0}) \ A^'. Notice that we have 

y 

dre{T, y, 0) = Kj {S,y) with S,y = --- 

T 

If 2/ > -S.J, then a^e(T, y, 0) < and T{y, 0) = 1. 

If 2/ < -C7, then ^^^ = ^J is a solution of dre{T,y,0) ^ Kj(^y) = and T{y,0) = -J^. This 
ends the proof of the lemma. D 
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Proof of Lemma \6.1(^ The proof proceeds in several steps. 



Ste p 1: c ontinuity. From Proposition 6.12 we already know that I'-'unetion ^ (^^((J* x J*)U A^*) 
and (6.191 holds true with 



where r = T{y,x) in the first line. Therefore, in order to prove that 'D'!^^^^^^^ G d^{Jj x J^), it 
is sufficient to prove that if {y,x) e {d{Jj x Ji)) \ A-?'^ = (({0} x J,) x (J, x {0})) \ A-'% and if 
{y'',x'^) e (J* X J*) \ A-'' is a sequence of points such that {y'^^x'^) — 7> {y,x), then we have with 



Tiy^x'^) 



£„ and — 



^. 



where {^y,^x) is given by (6.20). Let us recall that t*^ is characterized by the equation 



K, 



K, 



1-T* 



(6.37) 



(6.38) 



We will assume (up to extract a subsequence) th at r^ — >■ tq for some limit tq G [0, 1]. Because we 
have jx'^'p + ly'^p < B?, it is easy to deduce from (6.38), that there exists a constant C/j such that 



y 



<Cr 



(6.39) 



1 -rfe 
This can be proved by contradiction, distinguishing the cases tq — 0, tq — 1 and tq e (0, 1). Up 



to extract a subsequence, we can then pass to the limit in (6.38) and get 



K, i^y) = K, (£,) with ^y<0<^x 



(6.40) 



In the following cases, we now identify one of the two quantities ^y or ^j., the other one being 
determined by ( |6.40 ). 



Case 1: {y,x) e ({0} x J*) \ A-^'. From Lemma 
then 



6.15 



we know that tq = max f 0, 1 — ^ ) i s-nd 



and from (6.35), we get 



£, = max(a:,£+), Cy = (Kry^K^iC,)) 



^/u„ctio„(0,2:) = - (L.(?.) - io(0)) + io(0) 



e 



Case 2: (y, a:) G (J* x {0}) \ A^\ From Lemma 



(6.41) 



then 



6.15 



we know that tq = min ( 1 , — -^ ) , and 



and from (6.35), we get 



-e, = max(2/, -£-), £. = {K+)-\K,{^y)) 



y 



I?f,L„.„Jy, 0) = -f (L,(£,) - Lo(0)) + Lo(0) 



'junction Vy 7 ^/ ^ V J \^y / 



(6.42) 



Case 3: (y, x) G {(0, 0)} \ A-'*. This case only occurs if j G /q or i e /q. Moreover at least one of 

the two quantities — ^ and jz^ tends to zero. 
If £j, = 0, then Ki{£^,j.) — Kj{Q) and hence 



ey = o 



L,(0)>L,(0)=Lo(0) 
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If Cx = 0, then Kj{^y) = Ki{0) and hence 

e, -0 =^ Lj(0)>L,(0)=Lo(0) 

This imphes that 

^. = ^^+-0, Cy = ^i<0, if L,(0)=io(0)<L,(0), 
e. = C+>0, Cy-e^-O, if L,(0)>ij(0)-Lo(0), 
e. =e+ = 0, ^y = ^i=0, if L,(0)=ij(0) = Lo(0). 



By the uniqueness of the hmit, this finally shows that 2?^u„ction ^ ^^i^j ^ "^0 ^^^ (6.20) holds. 



Step 2: checking (6.35| and (6.22). From (6.41) and (6.42), we deduce (6.35) on ((J* x 



- TlJ' 



{0}) U ({0} X J*)) \ A^\ From I?^„,,,„„ = P^^^^^^ 

((j;x{o})u({o}x j;))nA^\ 



on A^% we deduce that (6.35) is also true on 



Then ( 6.22 ) follows from a simple computation for (y, x) ^ (0, 0) . This is still true for (y, x) — 0, 
because T^J^^tion i^ ^^ ■ This ends the proof of the lemma. D 

6.5 Study of "Dstraight 

The following lemma will be used below. Since it is elementary, its proof is omitted. 

Lemma 6.18 (Properties of I^gtraight)- ^°'' j ~ i ^ In, we have for (y,x) G Jj x J,; with 
(2/,a;)^(0,0) if j ^ i e In\ lo-' 

^straight (2^' 2;) - xd^Vll^-^^^.iy, x) - ydyVH^-^^^iy, x) = Lq{0) + K,{x - y) 

- Hj{-dyVil.^^{y,x)) 



d.V- 



iJ« 



i{y,x)^L[{x-y), dyVl'^^ {y,x) = -L'^{x-y) 



^ straight V 

6.6 Proof of Theorem 16.21 

We are now in position to prove Theorem |6.2[ We prove several lemmas successively. 
Lemma 6.19 (Properties of Pg*). For {y,x) £ Jj x Ji, we have 



Vi\y,x) 



Li{x - y) 

■^junction (.y' ^) 



if i= i & lo, 
if i + J, 



min(^j^„ction(y' x),L^{x-y)) if i = j E In \ lo- 



in particular Dq^ G ^'^{Jj ^ •^i) *" ^^^ fi'^'^^ ^^o cases. 

Lemma 6.20 (Singularities of the gradient of Dq). In the case i — j £ In \ Iq, we have 



Kiy^^) 



linear 



in a neighborhood of {d{Jj x Jj)) n A^*, 
Li(x — y) in a neighborhood of {d{Jj x Ji)) \ A^*; 



(6.43) 



moreover, in this case there exists a curve F^* such that T>q^ G C^{{Jj x Ji) \ (F-'* U {Yj, Xj})). 
This curve connects Yj — (—■^7,0) and Xi = {0,£,i) and is contained in (J* x J*) n A-'* 

The results of these two lemmas are illustrated in Figures [3] and |4] 
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A y ^ h 



j 7^' 



, ' 

y X 



yeJi 



Figure 3: i E Iq or j E Iq- trajectories never stay at the junction point 



n y^Jj 



kv^-h 



i^ti 



f*; 



y 




Figure 4: i, j E In\ Iq'- trajectories do stay at the junction point if (y. x) ^ (0, 0) 



Proof of Lemma 6.19 We only have to treat the case i ~ j ■ The convexity of Li imphes that for 

re (0,1): 

e(T,2/,x) = tL., y-'-j + (1-t)L,; (^37) - ^»(^-y)- 
Therefore for {y,x) £ J* x J* with j = i, we have 

^/mpiicitCy' 2;) == ^ mf^^ e(T, y, x) > L,{x ~ y) 

When i^je /q, we have 2?j'u„ctio„(y' ^) = ^hnpiicitl^' 2;), and then 

^j'unctionly.a;) > Li(a; - y) = I^straightC^'a;) = Vff\y,x) 
for (y,x) G J* X J* and then also for (y,x) G Jj x J^, by continuity of the functions. 
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Proof of Lemma \6.2C\ We first prove (6.43) and then describe the curve T 



-3,1- 



Proof of (6.43). Combining (6.131 and (6.17), we obtain 
Vl\y,x) = niin(2?j'^„^ti^„(?/,a;),I?^jVgi^j(2/,a;)) 

^ [^straight (2/' ^) = L^{x - y) for {y, x) e (J, X J,) \ A^-' 



(6.44) 



On the other hand, we have (a strictly convex function being above its tangent) for x ^ ^^ and 

y + ~q 

L^ix) > i.(e+) + {x - e+)iKe+) = xL'.i^t) + io(0) = 2?fLar(0,a:) 
L,{~y) > L,{^-) + i-y - ^-)L'^i^j) = -y^i^-) + Lo{0) = V^^^ijj, 0). 
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This shows that 

straight linear 



^^Uht > ^^Lar on (9(J,xJ,))nA^\ (6.45) 



We see that (6.441 and (6.45) imply (6.43) 



Description of {vf^^^ = V^^^^ } n A^"'. Notice that 

f 2?Lar(0,e+) = CtmCt) + io(0) = L,m = 2?^traight(0'^l')' 

1 ^fi;car(-c-,o) = e7i;(e7) + ^o(o) = L,(e-) = 2?^tWt(-e7.o). 

This means that the functions ^^f/ncar ^^d 2?g*j.g^j , ^ coincide at the two points Xi = (0, S,^) and 
Ij = (— ^~,0). Therefore we have 

^straight < ^unear ^n the Open interval ]X, , Y, [ 

because I'fincaj. is linear and T^H^^-y^^ is strictly convex as a function of y — a;. 

The function {y, x) ^ 2?straight(y' ^) " '^iLc^riv^ ^) being convex because X'j'4^^^(y, a;) is linear, we 
can consider the convex set 

K'' - {{y-.x) e J, X j„ PfLight(y,^) < 2?fLar(y,^)} ■ 

Then for i ^ j E In \ Iq, the set 

r^^ = {(y,x) e A^\ 2?Lar(2/,^) = ^ftraight(y>^)} 

is contained in the boundary of the convex set K^"^. More precisely, we have 

p* = {{dK^') n A^'^) c J* X J* 
which shows that F^Ms a curve and 

r^ = P''U{X,,Y,}. 

n 

Lemma 6.21 (The equations for Dq). For all i,j, and x,y where 2?g* is C^ : 

= -H,{{d,Vi'){y,x)) = -H,{{-dyVi'){y,x)). (6.46) 

Moreover for all x E Ji (with x ^ ^J' if j ^ i E In \ la) 

{VI' - xd.,V'^ - vdyVl'){Q, x) = LoiO) + if,(max(x, O) (6.47) 

and for all y E Jj (with y ^ — ^-^ ifj = iElN\ Iq) 

iV'^' - xd^V^' ~ ydyV'o'){y,0) ^ Lo{0) + K,{-max{y.-^-)). (6.48) 

We also have 

d^Vi\y,x) - L[{^,), dyV^\y,x) = -L'^{^y) (6.49) 

for all {y, x) E d{Jj x Jj) except for i = j E In \ Iq for which we exclude points {y, x) E {K,, Xi}. 
Moreover for j ~ i E Iq, we have 

Cx= ^y = x-y for all (y, a;) E 9( J, x J^) (6.50) 

and j — i E In \ Iq, we have 

iy = x~y, C.^x~y 

?, -^7, e. =e+ for iy,x)E{d{J,xJ,))nAJ^ ^''■^'^ 
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Proof. Using Proposition 6.12 for Vl^^^^twu^ Lemma 6.18 for I^gtrajghtj and (6.32) for Vl^^^^^ and 
the property (6.43), we geT 



{Vi'-xd,Vi'-ydyVi'){{),x) 



Lo(0)+ir,(max(a;,e+)) 

Lo(0)+if.(x) 

Lo(0) + X,(x) if x>i-^ 
Lo{0) if x<^^ 



which imphes (6.47). Similarly we get 



m'-xd,K~ydyK)iy,0) = 



Lo(0) + i^,(-max(2/,-C-)) 
Lo{0) + K,{-y) 

Lo{0)+K,{-y) if y>-i- 

Loio) if z/<-e7 



if i^j 

if i = j e Iq 

if i^ j e In\Io 



if i 7^ j 

if i^jelo 

if i^ j e In\Io 



which implies ( |6.48[) . Relations (6.46) and ( |6.49 ) follow both from Proposition 6.12 and Lemma 
6.18 Finally (6.50) and (6.51) follows from the previous results. This ends the proof of the 



lemma. 

We now can check the equations satisfied by Vq at the boundary. 
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Lemma 6.22 (Boundary properties of I?o)- ^t any point {y, x) e {0}x J^ with x^^^ifiE In\Iq, 
we have for any j € In 



{Vf;-xd.,Vf;~ydyVi'){y,x) = 



-ui&^kei«H-(-dyVf{y,x)) if N>2, 

~H^{-dyV^^{y,x)) if N = l. 



(6.52) 



Lemma 6.23 (Boundary properties of 2?o (continued)). At any point {y,x) £ Jj x {0} with 
y ^ -^7 */i ^ In\Io, we have for any i G In 



{Vl' - xd^V^Q - ydyVl'){y, a;) = - max H^. (d^Vl^'^iy, x)). 



(6.53) 



Proof of Lemma \ 6. 2S\ We first remark the general fact that 



On the one hand, from Lemma 6.19 we have for points [y^x) e {0} x Jj where Pq* is C^ 
-(2?S* - xd.^V'^^^ - ydyV^^){y,x) = Hk{-dyV',\y,x)) > H^ {-dyV^\y,x)) 



and this common quantity is independent on k. Therefore to conclude to (6.52 ) in the case N > 2, 
it is enough to show that there exists at least an index k such that 



Hki^dyV^,\y,x)) = H-{-dyVf{y,x)). 

Case 1: N > 2 and k ^ i. Then we have £,y < and then 

Hk{-dyV^^y,x)) = Hk{L'^{£,y)) = H^ [L'^Hy)) = H^i-dyV',\y,x)). 



(6.54) 



Therefore (6.54) holds true for k ^ i. 



Case 2: A^ = 1 and k = i = 1 E Iq. Then we have 

T^a{y,x) ^Vl^{y,x) ^ Li{x - y) 
and by Lemma 6.18[ we have for 



(Po" - ^d^T^l' - ydyVl'Ky, x) = -iJi {~dyVl\y, x) ) 



which is in particular true for y — Q. This shows (6.52) in the case iV = 1 
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Proof of Lemma \6.23\ From Lemma 6.19 we have for points {y,x) E Jj x {0} where Vq is C^ 
-iV^' ~ xd^Vi' - ydyV^'')iy,x) = Hk{d^,Vi\y,x)) > H^ {d^Vi\y , x)) 



and this common quantity is independent on k. Therefore to conclude to (6.53), it is enough to 
show that there exists at least an index k such that 



-ijfc 



Hk{d^Vl''{y,x))=H-{d^Vl''{y,x)) 



jk. 



Case 1: j € Jq. Then from Lemma |6.19[ we have with k ^ j 

dx'Dl^iy, x) = L'f,{S,x) with ix^x-y <Q. 



(6.55) 
(6.56) 



Therefore (6.53) holds true for k = j. 



Case 2: j £ In\Io- We distinguish sub-cases 
Subcase 2.1: y > -C~- From Lemma 



(6.53) for k = j. 



6.19 



Subcase 2.2: y < — ^- . Then we choose an index k £ Iq, and Lemma 

a,2?i^"(y,x) = 4(e.) with ex = 4+-o 



we still have (|6.56|) with k = j, which again implies 

implies that 



6.19 



which again implies (6.53) for such k (£ Iq. This ends the proof of the lemma. 
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We can now prove Theorem 6.2 



Proof of Theorem \6.S\ From Lemma 6.19 we know that 2?o has the regularity Cl except on certain 
curves F^* U {Yj.Xi} for j = i £ In\ Iq- So if {y, x) is a point of local Cl regularity of Vq, then 
we simply set 

ipo — Vq locally around {y,x). 

If (y, x) is a point where Vq is not Cl, then we have ^^{y, x) ~ 2?junction(y, x), and we can simply 
set 



^0 



-^junction ^^ ^ • 



The required equalities follow from Lemmas 6.19 6.22 6.23 Estimate (6.6) follows from the fact 
that 2?o is the minimum of ^junction G Cl{J'^) and of functions in C^{Jf) for some i. This ends 
the proof of the theorem. D 



A Appendix: Stability and Perron's method 

This section contains classical results from viscosity solutions, whose statements are adapted to 
the equation studied in the present paper. 



1.1 Stability results 



In view of Proposition 3.3 the following stability results are classical in the viscosity solution 
framework. See for instance [7]. 

Proposition A.l (Stability). Assume (AT) and let T > 0. 



• Consider a family of sub solutions (resp. super-solutions) {ua)aeA of (1.1) on Jt such that 
the u.s.c. (resp. l.s.c.) envelope u of 



sup Uo 

aeA 



(resp. inf u^ 



is finite everywhere. Then u is a sub-solution (resp. super- solution) of (1.1) on Jt 



• Consider a family of sub solutions (resp. super-solutions) iue)ee{o,i) of (1.1) on Jt such that 
the upper (resp. lower) relaxed semi-limit u is finite everywhere. Then u is a sub-solution 
(resp. super- solution) of (1.1) on Jt. 
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1.2 Perron's method 



In this subsection, we state the existence of a sohition of ( 1.1 1-( 1.2 ) which can be constructed by 
using Perron's method. This method is the classical way to get existence in a viscosity solution 
framework. 

Theorem A. 2 (Existence). Assume (AO)-(Al') and let T > 0. Then there exists an upper semi- 
continuous function u: [0,T) x J — )• M which is a viscosity solution of (1.1|-(1.2) on Jx and 
satisfies 

\u{t,x)-uo{x)\<Ct for (i, x) e [0, T) X J 
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